Unwinding in Hopfion vortex bunches 
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We investigate the behaviour of parallel Faddeev-Hopf vortices under energy minimization in a 
system with physically relevant, but unusual boundary conditions. The homotopy classification is 
no longer provided by the Hopf invariant, but rather by the set of integer homotopy invariants 
proposed by Pontrjagin. The nature of these invariants depends on the boundary conditions. A 
set of tightly wound parallel vortices of the usual Hopfion structure is observed to form a bunch of 
intertwined vortices or unwind completely, depending on the boundary conditions. 

PACS numbers: 11.27.+d, 05.45.Yv, ll.lO.Lm 



I. INTRODUCTION 



The standard model with knot soUtons, the Faddeev- 
Skyrme (FS) model, was proposed by L. Faddeev in 1975 
Since then, many analytical and numerical results 
have been obtained about this model 0, S, 0, H S 0, 
H, The localized solutions of the FS model are 

characterized by the Hopf charge, which can be defined 
if the field is constant at infinity: linv^oo ''^('r) = fioc', 
this allows one point compactification of the domain 
U{oo} ~ 5^ and the definition of Hopf charge using 
the homotopy group 773(8^). However, this is not the 
only possible boundary condition. Indeed, some physi- 
cal systems, like rotating superfluid '^He in the A phase 
[ll|,[l3| and certain insulating magnetic materials, called 
topological insulators can have different boundary 
conditions requiring a new topological analysis. There 
have already been experimental observations of a two di- 
mensional topological insulator in Bii_3;Sb:E, having the 
topological charge characterized by of Z2 x Z2 

In this paper we use large scale numerical optimisation 
routines to investigate the minimum energy configura- 
tions of the FS model in novel topological situations. The 
topological invariants that are relevant when the bound- 
ary conditions are no longer the usual limr^oo n{r) = ftcc 
are discussed in section[lll It turns out that new topolog- 
ical invariants, introduced by Pontrjagin are needed 
whenever some periodic boundary conditions are used. 
In section IIIII we describe the Faddeev-Skyrme model 
and the initial configurations. In Section [IV] we discuss 
the process leading to minimum energy configurations 
from initial configurations consisting of tightly wound 
and packed vortices parallel to the z-direction. We will 
demonstrate how the topological invariants of section lUl 
are conserved while allowing the unwinding of the Hop- 
fion vortices. The behaviour turns out to depend not only 
on the boundary conditions, but also on the size of the 
periodic cell as compared to the region occupied by the 
vortices. Intermediate states of the energy minimization 
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process are used to illustrate the processes that unwind 
the knot. Locally, these processes are quite similar to 
those seen in the context of closed Hopfions 



II. TOPOLOGY 

We study the FS model in situations where the phys- 
ical space cannot be compactified to but is periodic 
in some direction (s). In particular we consider the case 
when the domain can be identified as x T^ (periodic 
in the z-direction) or T"^ (periodic in all directions). In 
these cases the field n becomes a map n : x T^ — > 
or n : T'^ — > S^, respectively. Here we differenti- 
ate between and T^: If there is a special point x* 
such that f{x*) = f{—x*) (e.g. the point at infinity) 
we use while if the periodicity can be expressed by 
/(.t) = f{x + L), Vx then we use T^. 

Note that computationally, the difference between the 
three cases above is only that of the boundaries. If we 
want to model the S'^ boundary condition limr^oo n{r) = 
Hoo then in the computational lattice the fields at the 
edges of the box will be fixed to Hoo and no interaction 
can occur across the edges. Next suppose we use the same 
initial confinguration, but during minimization we only 
enforce the periodicity condition n(x™''^) = 7?(xf"") for 
some i. If the behavior of the system is such that the 
values of the field on the edges of the box do not change 
much the system will behave essentially as it behaves 
when the S"^ boundary condition is used. As we shall 
see, this is also refiected by the fact that in the T-^-case, 
the homotopy classification of the system is still given 
by the Hopf invariant even though the domain is not 
topologically S'^. The above systems behave differently 
only when there is strong interaction across one or more 
edges of the periodic box. 

The homotopy classification of the above cases is more 
complicated than the usual case of S'^ — > and was first 
proposed by Pontrjagin A more modern approach, 
restricted to closed, connected, oriented 3-manifolds is 
due to Auckly and Kapitanski We will further re- 

strict their result to maps from x T^ and T'^ to S^. The 
following geometrical descriptions are due to Kapitanski 
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Let n : X — > S^. The topological properties 
of r? follow the topological properties of its restrictions 
n|S^ and n\T^. If we denote these restrictions by / and 
g, respectively, we have / : — *■ and g : — > 
S^. The topological properties of / are described by the 
usual homotopy classification of maps ^ S^, i.e., by 
712(5^) and therefore the relevant homotopy invariant is 
S2 = deg / G Z. The topology of g is more subtle. It turns 
out that for each value of S2 there is another invariant, 
here called Si, which is defined by the map g. It can be 
shown that si G Z2S2 ■ 

As an example, consider the stereographic coordinate z 
of and the map n : x ^ for which f{z) = z'", 
where m G N. It can be shown [l^ that for such a map, 
deg / = m. Thus S2 is the number (with multiplicities) 
of zeros of /, i.e. the number of vortex cores. Now, for 
a simple vortex S2 = deg/ = 1 and therefore si G Z2. 
Thus, there are two homotopy classes, represented by 
ze^"^^ , where 6 G [0, 27r] is the coordinate of T^. The 
function g now describes how many full turns the is 
rotated as 9 goes from to 27r. The domain can be re- 
garded as a set of co-centric 2-spheres for different values 
of 6 with the inner- and outermost spheres identified. (If 
si — 0, g is the constant map 1— > and f{w, 9) — w, \/9 
and there is no rotation.) The full map n : x ^ 
represents an 1-vortex, where two distinct preimages of 
points on the target have linking number H , that we 
still call the Hopf invariant, but it is not a homotopy 
invariant in this case. 

Detailed description of how the higher numbers of ro- 
tations can in this case be unwound to either si = or 
si = I can be found in [l^, but we will not discuss this 
further. The case S2 = is rather special and is not 
needed here. 

Next, consider n : T'^ ^ S^. According to Pontrjagin 
[l5| . there are three primary homotopy invariants, ti, t2 
and <3. In order to determine the values of ti, consider 
n~^{p) G T^, the preimage of a regular value (under the 
map n) p G S^. By continuity of n, the preimage is 
a closed loop and therefore can be represented by a di- 
rected path 7 : 52 ^ Now, 7 G 7ri(T3) = Z x Z x Z 
and each of the invariants U belong to one of the Z of 
7ri(T'^). Visually, ti count the number of times 7 travels 
around the respective of the T'^ and can be determined 
as follows: For each ti, find all the points (with multiplic- 
ities) where 7 pierces a plane defined by Ci, numbering 
them with j G Z. For each of these points, let e(j) = -f-l 
if 7 is directed in the same direction as and e(j) = —1 
if the directions are opposite. (Note that by "piercing" 
we mean that the preimage has to go through the plane, 
not just touch it at a point.) Now, for each i, we define 

In the ft : ^ case there is also a secondary ho- 
motopy invariant, which further divides the maps with 
equal ti, t2, t^. For fixed values of ti, one defines 
t = gcd{ti,t2, ts) and for each t there is a secondary ho- 
motopy invariant, denoted ht. li t ^ 0, it can be shown 
that ht G Z2t [13, and thus for fixed ti and t > 



there are 2t different homotopy classes. If t = 0, this 
invariant is exactly the Hopf invariant and we denote it 
by ho. 

Two homotopic maps x ^ have the same in- 
variants Si and two homotopic maps T'^ — > have the 
same invariants ti and ht, but due to the secondary in- 
variants si and ht, maps with the same S2 or ti are not 
necessarily homotopic. In this work, we only deal with 
the primary invariants ti and S2. 

III. THE FADDEEV-SKYRME MODEL 

The explicit form of the Lagrangian of the Faddeev- 
Skyrme model can be written in terms of a unit 3— vector 

fi as 

C = C2d^,rfd^n + CAF^,,F^'\ (1) 
Ff,„ = ]^eabcn''d^,v!'d„rf, (2) 

where C2 and C4 are coupling constants. Choosing the 
usual metric (+, — , — , — ) yields the static energy density 

£ps = C2\\Vnf + c4Fjkf. (3) 

If one requires the field to have a finite energy in M'^, it 
is necessary to impose the boundary condition 

lim n(r) = ~ constant. (4) 

r — ^00 

However, if is the physical space is periodic in one or 
more dimensions, we have in the periodic direction(s) Xi 
n{xi) = n{xi+L), where L is the length of the periodicity. 
The space no longer compactifies to S'^ , but depending on 
the number of periodic dimensions, it becomes i) x T^ 
if one direction is periodic and the two others have the 
usual boundary condition Eq. (0]), ii) T^ x S^, in the case 
where two directions are periodic, and iii) T'^ if all direc- 
tions are periodic. The homotopy classifications of such 
fields were introduced in section [III We will not consider 
case ii) in this work. 

The initial configuration of a single vortex is the same 
as in [iOl. Using cylindrical coordinates p,9,z and two 
integers m, n, the form of the field ft is 

/ y/ l - f{p)^ cos{m9 + 2TTnz/L^)\ 
^ = - fipf sin(^^' + 27rnz/i,-) , (5) 

V f{p) I 

where Lz is the box length in the z-direction and / is 
just some profile function with /(O) = — 1 and f{p) — +1 
when p = 00, which is then cut into a finite size box. 

Case i) with initial states composed of a single straight 
tigthly wound vortex of ^ was studied in [20|. Those 
initial states correspond to S2 = 1 and values of the Hopf 
invariant larger than one. In principle they are homo- 
topic to configurations of lower Hopf invariants, namely 
H G {0, 1}, but there is an energy barrier which prevents 
the unwinding. 
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In this paper, the initial configurations are such that 
the preimages of (0,0,-1) are the cores of parallel vor- 
tices in the ^-direction; one such vortex is displayed in 
Figure [TJ The small bend visible in the figure was intro- 
duced in order to speed up the minimisation process and 
does not affect the topology or linking number of the sys- 
tem. Such an initial configuration has S2 = m in the case 
i) and ti = t2 — 0, = m in the case iii). In both cases, 
the linking number of preimages is mn. All of these are 
easy to verify visually. 




Figure 1: Single vortex in x or domain. The 
vortex core ns = — 1 and one point from the latitude 
ria = 0.1 are plotted as narrow blue and red tubes, 
respectively. The wireframe indicates the boundary of 
the box. 



Vortex bunches are constructed by putting several 
identical vortices into the computational box, as shown 
in Figure [Sal For the multi- vortex configurations inves- 
tigated here, we have used the same values n = 6 and 
771 = 1 for each vortex. The value of n may seem rather 
high, but we have learned in [131 that for low values the 
the vortex does not bend much and therefore would not 
properly interact with the neighboring vortices in the cur- 
rent situation. 

In this work, we use the same programs as in [3, 0, [3, 
[20l | because the only difference in the computations are 
the boundary conditions. For numerical computations, 
this change of boundary conditions is rather minimal, 
since the MPI library directly supports both periodic and 
fixed of boundary conditions defined independently for 
each dimension. 



IV. RESULTS 

We now describe how an initial configuration of a 
bunch of Hopfion vortices built from Eq. ^ continu- 
ously deforms to a minimum energy configuration while 
conserving the homotopy invariants Si, U. We will also 
describe, how, depending on the details of the initial con- 
figuration, the linking number or the Hopfion (the Hopf 
invariant) is sometimes conserved and sometimes not. 

All the computations have used fully periodic bound- 
ary conditions, but in different computational lattices. 



A. 2x2 vortices in a large box 

Let us first consider an initial configuration consist- 
ing of 2 X 2 vortices parallel to the z axis, put into a 
very large box. It is expected that there would be neg- 
ligible interaction across the edges of the box, i.e. the 
system would behave as it would in x T^. In the ini- 
tial configuration, it was necessary to pack the vortices 
into a tight formation with enough twisting {n = 6) to 
ensure that they interact, but not across the boundary. 
Technically, this was accomplished by building the ini- 
tial configuration from 16 narrow boxes: the four central 
ones contained a single vortex each, while the remaining 
12 pieces were filled with the vacuum, n = (0,0,1). In 
contrast to we do not fix the field value at the x, y 
boundary but only require periodicity, this allows us to 
see whether the behavior of the vortex bunch is the same 
in x T^ and T^ x T^, where the T^ is the large box in 
the xy-p\ane. Indeed, the values on the xy-boundaries 
do not change much during the minimisation process. 

The initial configuration is shown in Fig. [Jal Since 
= (0,0,1), as seen from Eq. ([5]), the vortex core is 
at 77 = (0, 0, —1). It is easy to count the linking number 
of preimages, which equals to 24 (recall that there is a 
^ factor involved - there are 48 crossings) and also the 
homotopy invariants U: ti and t2 are zero (no winding 
the in x- and y-directions) and — A (four cores pierce 
the xy-p\ane). 

Energy minimisation is applied to the initial config- 
uration using a gradient based algorithm 0|. Initially, 
the interaction between vortices is negligible, but after 
around 40000 iterations, certain parts of different vor- 
tices are close enough to produce evolution which differs 
from that of a single vortex. At 50000 iterations, these 
parts touch and the corresponding preimages reconnect, 
joining two vortices together. Fig. HH shows the situation 
before and [2c] just after the reconnection. The wireframe 
box showing the initial locations of the vortex cores is 
also displayed for reference. The reconnection process is 
allowed because mathematically, the red curve is a single, 
albeit multiply connected, preimage, and the elementary 
process is exactly the same as is described in [l3| for a 
single vortex. Such deformation processes repeat several 
times in different parts of the lattice until after 300000 
iterations, the initially separate vortices have formed a 
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tangled bunch where the cores (and preimages of a point 
where = 0) wind around each other (see Fig. \2d\i . 
Visual inspection of the intermediate configurations con- 
firms that, as expected, there is no interaction across the 
xy-boundary. Note that there are six cores (and preim- 
ages) piercing the the box at the top and bottom edges, 
but the direction of one of them is opposite to the others, 
giving ^3 = 4 as required. 

The final, minimum energy configuration, has the same 
linking number of preimages and integers ti as the initial 
state. This is despite the fact that the linking number is 
no longer guaranteed to be conserved. Its conservation 
is solely due to the boundaries of the box being very far 
from the vortices. Indeed, in the next part, we will see 
that nonconservation is possible. 



B. Vortices in a small box 

In this case the same initial configuration is put into 
a box half the previous size, with the expectation that 
there now would be interaction across the edges of the 
box. Technically, this was accomplished by taking the 
same four vortices as above but omitting the 12 vacuum 
pieces. Again, the x, y-directions are periodic but now 
the vortices are equally spaced with no change in the 
spacing at the boundary of the computational box. 

The beginning stages of the energy minimization pro- 
cess are identical to that of the large box, but when the 
vortices start to feel each other the configuration de- 
velopes in a different manner. Instead of reconnecting 
to join two vortices together, the reconnected preimages 
tend to develop into straight filaments extending through 
the whole periodic lattice. Later, some of the filaments 
again reconnect and form loops, but since these loops 
have trivial linking with all other preimages, they can be 
deformed into nothing. 

Thi s is illustrated schematically in Fig. [3l where ini- 
tially (a) , processes familiar from [l^, [lol have already 
led to each of the four preimages to split into a relatively 
straight part parallel to z-axis and a loop around it. The 
preimage loops extend, eventually touching each other, 
and then reconnect along the dashed curves |(b)[ Upon 
further energy minimisation, the filaments will in turn re- 
connect along the dashed curves to form loops |(c)[ The 
resulting loops have trivial linking with other preimages 
and are threfore free to collapse and vanish [(d)| The ver- 
tical preimages stay relatively unchanged through this 
process. This is the general mechanism of Hopfion un- 
winding. 

Indeed, the energy minimization process shows how 
this process of filament and loop formation, followed by 
the loops shrinking and vanishing, is repeated a number 
of times in different parts of the computational box so 
that the final linking number of preimages is zero. The 
only preimages left in the system are preimages piercing 
the ccy-plane. At the same time, the energy approaches 
the energy of one or more straight unwound vortices. 



In order to calculate the invariants ti during the min- 
imisation process, it is necessary to keep track of the 
directions of the preimages, just as must be done for the 
computation of the linking number. Now, initially we 
have ti = t2 = and = 4, so whenever a preimage 
reaches a boundary of the computational box and pierces 
it, there must always be a corresponding piercing of the 
same boundary but in the opposite direction, so that ti 
are unchanged. Indeed, this has been observed. 

In order to further investigate the processes involved 
in the unwinding of the Hopfion, we simulated a system 
where the initial condition is a single vortex. Again, we 
use m = 1 and n = 6, but instead of packing several 
vortices together, we take just one vortex ^ and add 
no vacuum to the system. The boundary conditions are 
again fully periodic, but now the length of the period 
in xy-directions is half of that in IIV Bl Any splitting 
and recombination of preimages is now expected to occur 
across the lattice boundary. This is indeed observed with 
no new processes present in the system. 



V. CONCLUSIONS 

We have studied how the Faddeev-Skyrme model be- 
haves when the domain of the model is taken to be 

X T^ and T^ instead of the usual S"^. The new do- 
mains require a new homotopy classification which is due 
to [H, Il3, For X T^ there is a primary homotopy 
invariant S2 related to the part of the field configura- 
tion and when S2 ^ 0, there is also a secondary invariant 
si related to the T^ part. We discuss the invariant S2 
only and it is seen to be conserved in the energy minimi- 
sation process. For T'^, there are three primary homo- 
topy invariants, ti, t2, ^3 related to the winding numbers 
around the three tori T^ of the T^. There is a secondary 
invariant in this case also. When at least one ti 0, 
one defines t = gcd(ii, ^s) and the secondary invariant 
takes values in Z2t- If all ti = 0, the secondary invari- 
ant becomes the Hopf invariant. We have discussed the 
primary invariants only. 

We have constructed a set of initial configurations con- 
sisting of a number of twisted vortices where energy min- 
imisation is known to lead to knotted Hopfions in the 

X T^ case [23|- It is observed in this case, that in the 
continuous deformation driven by energy minimization 
the conservation of si and S2 does not allow for the un- 
winding of the knot and leads to an intertwined vortex 
bunch. In the T"^ case neither the conservation of ti nor 
energy minimisation requires that the knotted structure 
remains intact and indeed the Hopfion unwinds. The 
conservation of the homotopy invariants is confirmed by 
a visual inspection in all cases. 

The results have been compared with previously known 
deformation processes of preimage splitting and recon- 
necting and found to follow the same deformation "rules" 
as previously [iB, US]- No new deformation processes 
were observed. 



(a) Initial state of 4 vortices with n = 6 each 



(b) The preimages are about to touch at (50000 
iterations) 





(c) The reconnected preimages have joined to 
vortices together (60000 iterations) 



(d) Finally, (after 300000 iterations) the whole set of 
vortices is bundled into a single bunch by individual 
preimages wrapping around all the cores. 



Figure 2: Relaxation of an initial configuration of four vortices in x domain. The colours are as in Fig. (TJ but 
the wireframe now shows the size of the lattice in z-direction and the exact locations of the vortex cores (the vertical 
wires) . The computational box is larger and the thicknesses of the tubes are arbitrary. 













c 































(a) 

Figure 3: Schematic illustration of the formation of filaments and contractible loops. The solid and dashed lines 
represent the preimages and the dotted parallelepiped is the outer boundary of the periodically repeating region 

containing the four unit cells of the initial vortices. 



The results presented here have possible experimental can have topologically nontrivial properties (l3l |. 
relevance to the observations of vortices in a rotating 
superfluid ^He in the A phase Another physically 

relevant situation where maps fron T'^ arise are 

insulating magnetic materials in three dimensions, which 
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